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Abstract 

Considering the unphysical result obtained in the calculation of the free-energy cost 
for twisting the boundary conditions in a spin glass, we trace it to the negative 
multiplicities associated with the Parisi replica-symmetry breaking (RSB). We point 
out that a distinct RSB, that keeps positive multiplicities, was proposed long ago, 
in the spirit of an ultra-long time dynamical approach due to Sompolinsky. For 
an homogeneous bulk system, both RSB schemes are known to yield identical free 
energies and observables. However, using the dynamical RSB, we have recalculated 
the twist free energy at the mean-field level. The free-energy cost of this twist is, as 
expected, positive in that scheme, as it should. 
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1 Introduction 



Lately a rather strange result was uncovered. As is well-known for an 0{N) fer- 
romagnet, the breaking of the continuous symmetry and the associated Goldstone 
modes, tax the forcing of a twist in the orientation of the magnetization (between 
the z = and z = L boundaries) by an amount of free energy proportional to 
L'^~^fT . Consequently the lower critical dimensions above which symmetry break- 
ing occurs at a non-zero temperature is exactly given by dc = 2. For a spin glass, 
the broken continuous symmetry group is the reparametrization (or gauge) group. 
The analog of a rotation twist, is now a (small) gauge twist between the z = and 
z = L boundaries. In contradistinction with the 0{N) case, through a long and 
arduous calculation making use of Parisi's replica symmetry breaking (RSB) [2 on 
Parisi's truncated Hamiltonian, we obtained a twist free energy cost [3j proportional 
to —L'^~'^^^ (i.e. with a negative coefficient) |3] ; the exponent rj is the usual order 
parameter anomalous dimension, computed there to one loop. The implication for 
the lower critical dimension, namely (ic = 2 — ~ 2.5, was indeed in agreement 
with previous estimates 01^]. However the negative sign of the coefficient, i.e. a 
gain in free energy under twist, was very puzzling. Taking this at face value could 
point to the instability of a solution with a space inhomogeneous order parameter. 
If one considers for instance an Ising antiferromagnet, one may obtain a lower free 
energy with twisted, i.e. antiparallel boundary conditions. But, given that Parisis's 
space-homogeneous solution is now proven to give the exact free energy [71 [H] , this 
seems unlikely. 

A way out of this puzzle seems to be the following. Parisi's solution is (semi) 
stable, all eigenvalues of the associated Hessian being non-negative when the number 
of steps R of RSB goes to infinity 9j. However the multiplicities of those eigenvalues 
are all negative. This means that the saddle-point where the free-energy is calculated 
has the characteristics of a maximum ; hence small excursions away from it will yield 
an unphysical negative free-energy cost. This situation, in turn, is due to the fact 
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that, in Parisi's RSB, the natural ordering of the box sizes Pu (where u is the discrete 
overlap) is reversed, i.e. one works with 

Pr+1 = I > Pr > ■ ■ ■ > Pu > ■ ■ ■ > Pi > Po = n. 

In this article we would like to reconsider this question of the free-energy cost under 
a twist of boundary conditions, in the light of a distinct RSB scheme proposed long 
ago^ni- This alternative scheme has the following characteristics : 

• The box order is not reversed, hence one has multiplicities remaining 
positive: 

PO > Pi > ■ ■ ■ > J9„ > ■ ■ -pi? > Pr+i = 1 

• Whereas pu in Parisi's approach is a parameter without any physical 
conjugate field, it is replaced here by a susceptibility derivative A„, 
a physical quantity associated with a dynamical-like approach to the 
spin-glass problem (whose history is sketched below). In fact it is a 
better candidate for an order parameter, since it vanishes above the 
Almeida- Touless ^T] line and is non-zero below. 

• At the saddle-point it gives exactly the same value as Parisi's for all 
observables. It yields also the same eigenvalues for the Hessian at the 
saddle point |12j, but with different individual multiplicities. 

• The cost of a twist becomes now positive, at least to lowest order ( for 
the "kinetic" part of the free energy, i.e. the part which depends of the 
spatial variation of the order parameter, since the potential energy is 
invariant under the twist). The one- loop contribution will undoubtedly 
take some time to sort out. 

To replace in context the choice made here, a short historical reminder seems appro- 
priate. It all started with Sommers ^H] proposing a new solution with non-negative 
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entropy to the Sherrington-Kirkpatrick [TT model, a solution possessing a linear re- 
sponse anomaly. Enlarging on the RSB proposal of Blandin et al. soon after, 
De Dominicis and Garel fTB|, Bray and Moore |^ linked the non-negativity of the 
entropy to the infinite limit of what has been called above the box size po, a limit that 
restituted jT7| the Sommers solution. Substituting opaque replicas by dynamics [IS] 
led to some physical insight and Sompolinsky jTH] proposed to describe spin glasses in 
terms of the spin correlation function q = {era) and a spin susceptibility, or linear re- 
sponse X = (o'er) (o" being coupled to a magnetic field). The linear response anomaly 
A is closely related to x- the long time limit (when the initial time in Langevin's 
equations is sent to minus infinity) he described the sytem through an infinite set 
of relaxation times r„, with u = 1,2, - ■■ R and R —>■ oo. When u remains discrete , 
his g(r„) is what we call g„ and his A{t^^), our (with — A^ = A^ — A^+i). The 
approach used in the present article is thus a replica reformulation of Sompolinsky's 
dynamics as in JU]- At this point one may ask why not retain dynamics and drop 
replicas altogether. The answer is that 

• although some results have been obtained beyond mean-field through 
dynamics by Sompolinsky and Zippelius [201 HI] , if is more difficult to 
work with four time variables than with four replicas. For instance, the 
spectrum of masses (the eigenvalues of the Hessian) are easily obtained 
with replicas, but not yet fully sorted out within the dynamics. 

• Besides, one is interested in understanding how to obtain a physical 
answer for the twist free-energy. 

In the following we first compute in some details the contribution to the free energy 
cost of a twist in the parametrization, using Parisi's RSB, as already sketched before 
[H]. Then we perform the same calculation for the RSB introduced by De Dominicis, 
Garel and Orland ^| (DGO), that we might call also "dynamical-like" for lack of a 
better name, getting then the opposite sign. Finally we show that this result may also 
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be understood from the fact that there is a plateau contribution in Parisis's solution, 
where the overlap fuction remains constant and equal to its Edwards- Anderson value, 
and does not fluctuate. In the dynamical-like approach there is no such plateau. 
Hence when considering spatial variations the two approaches give distinct reults, 
whereas they lead to identical results for the bulk mean-field problem. 

2 Kinetic free energy a la Parisi 

In terms oi n x n matrices qab {qua — 0) the free energy, in Parisi's truncated model, 
reads for a spatially homogeneous order parameter 

nF(^)/L^ = - Y^C-ql + ^q^,) + | ^ g„,g,egc.. (1) 

ab abc 

The contribution of the kinetic part of the free energy Fk (i.e. the part which comes 
form a non-spatially homogeneous order parameter) is given by 

FP = ^ E E(?af (0 - ^ab{^ + (2) 
1=0 a,b 

in which we assume twisted boundary conditions in the z-direction (and periodic 
boundary conditions in the remaining (d-1) dimensions) ; space in the z-direction is 
kept here discrete with < i < L. Using Parisi's overlap function q^ , u — 0,1, ■ ■ ■ , R 
(with qaa = qR+i — 0), and the associated box size Pu (with po = n,pR+i = 1), one 
obtains 

I R+l 

-E^ab(^) = HiPuii) - Pu+i{i))ql{i) (3) 

^ ab u=0 

One needs also 

1 ^+^1 1 

- E (lab{i)qab{i + i) = E [9 +Pu{'i + 1)) - 7;iPu+i{i) +Pu+i{i + l))]g«(i)g«(i + 1) 

ab u=0 ^ ^ 

(4) 

This gives then 

j^d-l L-l R+l 

fP = — — E E[(^^«(0~^^"+i(0)?«(0~(^^"(^+i)~^^«+i(^+i))^«(^+i)][^«(0~?«(^+i)] 

^ i=0 u=0 

(5) 
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In the R oo continuum limit , we have 

Puii) - Pu+iii) -p{u; i)du (6) 

and finally 

Fr = du—{p{u]z)q{u]z))—q{u]z). (7) 

At the bulk saddle-point one has 

w 

= (8) 

For twisted boundary conditions, to lowest order in the twist h{u) ( /i -C 1 and 
limited to a small support < u < x), 

111 7 711 

q{u- z) = ^^{u + jKu)) + 0{h') = ^^p{u; z) + 0(/.2) (9) 
Substituting (0) into Q we finally obtain 



Jo ^""^ 

that is a negative cost in free energy for the introduction of a twist. 

3 Kinetic free energy a la DGO 

Whereas the zeroth step in Parisi approach is the {R = 0) matrix qab = q for a ^ b 
(with qaa = 0), here the zeroth step is the so-called Sommers starting point (formally 
identical to the -R = 1 Parisi's RSB). Namely, one divides the n x n matrix Qab into 
n/po blocks g^^/j (of size po xpo) along the diagonal, and ^(^ — 1) off-diagonal blocks 
ra,i3 (of the same size po x po) . Then one does R steps of RSB on both matrices qa,p 
and r^,/? (whereas in Parisi's approach those steps would only invlove the diagonal 
blocks qct,!^)- One thus needs a double labelling for each element of the initial matrix 
Qab in order to specify for a (and b) the po x Vo block matrix under consideration, and 
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the element in this block matrix . One thus writes a = {a, x) with x = 1, 2, ■ ■ ■ , n/po 
and a = 1, 2, ■ ■ ■ ,po 

Qab = Qa,x;/3,y = <la,f3^x,y + a,/3(l — ^x,y)- (H) 

The RSB steps apply now to the matrices qa^p ( with „ = 0) generating go, Q'l, • " " ; iRi iR+i = 
0, and to the matrices giving rise to tq, ri, ■ ■ ■ , r^j+i. The successive steps of RSB 
involve the box sizes po,pi, ■ ■ ■ ,Pr {pr+i = !)• In the DGO scheme the box sizes are 
made to go to infinity in the prescribed natural order 

Po > Pi > ■ ■ ■ > Pi? > PR+i = 1 (12) 

with at the same time letting the matrix elements of g — r go to zero in such a way 
that 

- Ao = po(go - ro) (13) 

remains finite and in general 

- A„ = Pu[{qu - ru) - (p«-i - ru-i)] ~ Pu{qu - ^n) (14) 
Here can be shown to be the (discrete) susceptibility derivative 

- A„ = A, - A„+i (15) 

(—Ay is a positive quantity). 

As a result the Parisi free energy functional F'^^^jp^; g„} is replaced by a dynamical- 
like free energy F*^^){A„; g„} with the associated stationarity conditions determining 
qu, A„ and their relationship. 

This alternative formulation gives in the R oo continuum and at the saddle- 
point lesults that are identical to Parisi's. Besides it also leads to the same eigenvalues 
of the Hessian at the saddle-point [T2j. 

After this brief reminder we now proceed to derive the kinetic contribution : 

4""^ = ^EEiQabi^ - QaS + 1)]^ (16) 
^'"^ i a,b 
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We have : 

^T.Qlb = ^ E - ra,f3)Scc,y + ra,f3? (17) 

''' ab "' a,l3 x,y 

Summing upon x, y one obtains 

lT.Qlb = ^ E{r-[(9",/5 - ^o.,^)' + 2^",/3(?-,/3 - V/?)] + (r-)'</3} (18) 

^ ab ^ a,f3 PO PO 

After R steps of RSB one obtains 

^T.Qlb = ^%o{Y.(p^-p^+^)[(il-''l) + -r'^n] 

^ ab ^PO u=0 PO 

= j:Pu[{ql-rl)-{ql^,-rl_,)]-rl^, (19) 

n=0 

Since p„ — >• oo, — >• 0, — r„ ~ ^ for m = 1, ■ ■ ■ , i? in the hmit (fT^ the equation 
(Hn)) reduces to 

' EQa6 = -2EgnA„-4+, (20) 
It 

The analog of Q is now 



^ ab ii=0 



r E Qlb{^ = -2 E ?.(^) A„(^) - 4+1- (21) 

Likewise (H)) becomes 
1 

E'5afe(«)<5afe(« + 1) = E^'«['?«(*)^«(^ + 1) ~ ru{i)ru{i + 1)] - (22) 

The box sizes going to infinity do not carry a space index any more and we end up 
with 

1 ^ . 

-Y.Qab{i)Qab{i + l) = E[^-(^)?«(* + l) + ^«(* + l)^«(^)] -4+1 (23) 

^ ab u=0 

The resuh for the kinetic part of the mean field free energy functional is thus 
= ^J:J:iQab{^)-Qabi^ + l)f 

j=0 ab 
j^d-l L R 

J2 [qu{i)K{i) + qu{i + + 1) - qu{i)K{i + 1) - qu{i + 1) A„(i)] 



2 

^ i=0 u=0 



' E E[(?-.(0 - qu{^ + i))(A„(0 - A„(. + 1))] (24) 



2 

^ i=0 u=0 



and in the double continuum limit, in which both space is continuous and R, the 
number of steps of RSB, goes to infinity 



At the saddle-point one has [101 CH] 

- A(n) = (26) 

w 

which is the " anomalous fluctuation-dissipation relationship" for the spin-glass. Again, 
to lowest order in the twist h(t), one can just keep, as in 

111 7 

q[u-z) = ^^[u+-h{u)]+0{h') (27) 



) = 1(^)2L'^-2 r duh\u) (2^ 



from which one obtains, to lowest order, the free energy for that twist 

4^2g' Jo 
so that 

fP = -2fP (29) 

The signs are opposite : the free energy cost under twist is positive for the dynamical- 
like RSB. 

4 A simple calculation at the saddle-point 

Given the previous difference between the two RSB schemes when one considers the 
spatial variation of the order parameter, one may examine in the light of the previous 
calculation why the two free energies happen to coincide for a spatially uniform order 
parameter. 

For instance let us consider the bulk contribution to the free energy per unit 
volume which is quadratic in the order parameter : 

ab 
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In Parisi's scheme one finds 

R-l 



^ u=0 

duq^iu) + (xi — l)g^(xi)] 
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= i(^)'[-y + (^i-i)^i] (31) 

Here the second term in the bracket comes from the " plateau " sector of q{u) 
{xi < u < 1) , at which it remains fixed to its Edwards- Anderson value 

w 

= -xi (32) 

The mean field stationnarity conditions determine xi in terms of the external param- 
eters (temperature and coupling constants). 
For the dynamical-like approach one finds 

^ u=0 



[2 / duq'iu)^qiu)-q'ix,)] 

4 Ju=Q W 

V2g' ^ 3 ^ ' 

The two results (jSD) and (jH^ at the end coincide as announced, however the origin 
of the various terms are different. Indeed in ()31|) part of the answer comes from 
the plateau where fluctuations are not allowed under twist. Thus when considering 
fluctuations the plateau part in ()31|) will not contribute and one will obtain 

(/(^))//.c. = -2{fP)fiuct. (34) 

as we found in (jSHI)- 

Similar calculations for the parts of the free energy which are cubic and quartic in 
the order parameter would show a similar conspiracy to make the bulk contributions 
identical, but yield different results when the order parameter varies in space. 
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5 Conclusion 



We have thus demonstrated that two approaches, using different RSB schemes, that 
provide identical results for the bulk mean-field energy, lead to distinct (and opposite) 
answers when one enforces spatial variations of the order parameter. The negative 
multiplicities occuring in Parisi's scheme were responsible for a decrease of the free 
energy under twist, whereas the dynamical-like DGO scheme does lead to an increase 
under twist. 

This has been established only at the mean-field level. It would be conforting to 
extend this to one- loop calculations at least, as we have done in a previous work j3] 
for Parisi's scheme, but the algebra, although well-defined, is quite elaborate. 

Given that the two RSB schemes yield the same bulk mean-field free energy, but 
are already different under twist, it becomes of great interest to compare the effect 
of higher loops on the bulk free-energy. 
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